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I ■■ Abstract 

I , Within framework of the fi from u Supersymmetric Standard Model (/ii/SSM), 

I exotic singlet right-handed neutrino superfields induce new sources for lepton- 

r-| I flavor violation. In this work, we investigate some lepton-flavor violating 

O-i' processes in detail in the /iz/SSM. The numerical results indicate that the 

Q^- branching ratios for lepton-flavor violating processes fi -^ 67, t ^ fi'y and 

/i — )■ 3e can reach 10~^^ when tan/3 is large enough, which can be detected 

in near future. We also discuss the constraint on the relevant parameter 

^ I space of the model from the muon anomalous magnetic dipole moment. In 

^ • addition, from the scalars for the /zi^SSM we strictly separate the Goldstone 

t;;;-j- ■ bosons, which disappear in the physical gauge. 
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.5^ ■ 1. Introduction 

X: 

5^ ■ It is obviously evidence of new physics beyond the Standard Model (SM) 

that if we observe lepton-flavor violating (LFV) processes in future experi- 
ments, because the lepton-flavor number is conserved in the Standard Model. 
In supersymmetric (SUSY) extensions of the SM, the R-parity of a particle 
is defined as i? = (—1) -^+3^+25' j^j^j g^j^^ q^qj^ \^q violated if either the baryon 
number [B) or lepton number (L) is not conserved J2|, |3| , where S denotes 
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the spin of concerned component field. Note that R = +1 for particles and 
— 1 for superparticles. 

Differing from the models in Refs.J2|, |3|, the authors of Ref. ^] propose 
a supersymmetric extension of the SM named as the "/x from u Supersym- 
metric Standard Model" (/xz/SSM), which solves the u problem |5| of the 
Minimal SupersjTiimetric Standard Model (MSSM) [6] through the lepton 
number and R-parity breaking couplings between the right-handed neutrino 
superfields and the Higgses eah^i^'^H'^H^ in the superpotential. The effective 
/i term eahl^'H'^H^ is generated spontaneously through right-handed sneutrino 
vacuum expectation values (VEVs), /x = Aj (i>f), as the electroweak symme- 
try is broken (EWSB). Note that a popular model is the so-called Bilinear 
R-parity Violation (BRpV) model |3| , where the BRpV terms eab^tH^LI are 
added to the MSSM. The effective BRpV terms are generated spontaneously 
through the R-parity conserved terms eabY^^^P'^H^L'l in the superpotential of 
the /iz/SSM, and Si = Yy^^ {^j)y ^s EWSB. So largely differing from the other 
models |2|, |3|, the /iz/SSM introduces three exotic right-handed sneutrinos z>f, 
and once EWSB the right-handed sneutrinos give nonzero VEVs. In addi- 
tion, the nonzero VEVs of right-handed sneutrinos induce new sources for 
lepton-flavor violation. In this work, we analyze the constraints on param- 
eter space of this model from the experimental observations on some LEV 
processes and muon anomalous magnetic dipole moment (MDM). 

If the left-handed scalar neutrinos acquire nonzero vacuum expectation 
values when the electroweak symmetry is broken , the tiny neutrino masses 
are aroused [TJ to account for the experimental data on neutrino oscillations 
[Sl, y, llO| . Three flavor neutrinos l'e,^J.,T are mixed into three massive neutrinos 
^1,2,3 during their flight, and the mixings are described by the Pontecorvo- 



Maki-Nakagawa-Sakata unitary matrix Upj^^^^ [Hj . The experimental obser- 
vations of the parameters in Up^^j^^ for the normal mass hierarchy 12| show 
that |13|] 

Am^i = 7.58lS-^^ X lO^^eV^ , Am^2 = '^■'^^-Vm x 10"^eV^ , 
sin^ e,, = 0.306+°:°}^, sin^ 9,, = 0.42+°;0«, sin^ ^13 = 0.021+^^8'- (1) 

Note that the Daya Bay Reactor Neutrino Experiment has measured a nonzero 
value for the neutrino mixing angle ^13 with a signiflcance of 5.2 standard 



deviations recently [lj|. Differing from the BRpV model, where one neutrino 



mass is generated at tree level and the other two at one loop [15|, the /iz/SSM 



can generate three neutrino masses at the tree level through the mixing with 



the neutralinos including three right-handed neutrinos [16|, ll7[. Here, we 
use the neutrino experimental data presented in Eq.([T]) to restrain the in- 
put parameters in the model. Then, we analyze the branching ratios for the 
various LFV processes: /i — )■ 67, r — )■ /i7, /i — >• 3e, etc., and the corrections 
to the anomalous magnetic dipole moment of the muon a^ in the /iz/SSM. 
The numerical results indicate that the new physics contributes large correc- 
tions to the branching ratios of the mentioned LFV processes and a^ in some 
parameter space of the model. 

The outline of the paper is as follow. In section |2l we present the in- 
gredients of the /xz/SSM by introducing its superpotential and the general 
soft SUSY-breaking terms, in particular we strictly separate the unphysical 
Goldstone bosons from the scalars. In section El we analyze the decay width 
of those interested rare LFV processes, and present the SUSY contribution 
to muon MDM in section HI The numerical analysis is given in section |5l 
and the conclusions are summarized in section [61 The tedious formulae are 
collected in Appendices. 

2. The /xi^SSM 

Besides the superfields of the MSSM, the /zz/SSM introduces three exotic 
gauge singlet neutrino superfields uf. The corresponding superpotential of 
the /iz/SSM is given as [^ 

- e^.Kv'rHlHl + \K,,,OtO'^^Ol , (2) 

where Hj = [hI H') , H^ = (i^+, H^) , Qj = («„ d^ , Lj = (i>„ e.) (the 

index T denotes the transposition) are SU{2) doublet superfields, and rf^, 
ifj and e^ represent the singlet down-type quark, up-type quark and lepton 
superfields, respectively. In addition, Yu^cL,e,vi A and n are dimensionless ma- 
trices, a vector and a totally symmetric tensor, a, h are SU(2) indices with 
antisymmetric tensor ei2 = —621 = 1, and i,j = 1, 2, 3. The summation 
convention is implied on repeated indices. 

In the superpotential, the first three terms are almost the same as the 
MSSM. Next two terms can generate the effective bihnear terms eab^iH^L"^, 
€abfJ'H°[Hl[, and Ei = Y^-. {^f}-, fJ- = Aj(z/f), once the electroweak symmetry 



is broken. The last term can generate the effective Majorana masses for 
neutrinos at the electroweak scale. And the last two terms explicitly violate 
lepton number and R-parity. 

The general soft SUSY-breaking terms in the /xz/SSM are given by 
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M3A3A3 + M2A2A2 + MiAiAi + H.c, 



(3) 



Here, the front two lines contain squared-mass terms of squarks, sleptons and 
Higgses. The next two lines consist of the trilinear scalar couplings. In the 
last line, M3, Mi and M\ denote Majorana masses corresponding to ^[/(S), 
SU{2) and U{1) gauginos A3, A2 and Ai, respectively. In addition to the 
terms from Cgoft, the tree- level scalar potential receives the usual D and F 
term contributions j^. 

When the electroweak symmetry is spontaneously broken (EWSB), the 
neutral scalars develop in general the vacuum expectation values (VEVs): 



Thus one can define neutral scalars as usual 
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For simplicity we will assume that all parameters in the potential are 
real in the following. After EWSB, the scalars mass matrices Mg, Mp and 
M|± are given in Appendix B The CP-odd neutral scalars mass matrix Mp 
contains a massless unphysical Goldstone boson G°, which can be written as 
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with an 8 X 8 unitary matrix Z^ 
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where v^^ = ^Jv\ + v\ and v^-^ = \/v\ + v\ + Vy^Vy^. Making use of the 
minimization conditions of the tree-level neutral scalar potential, which are 
given in 



Appendix A[ we have 



{ZlMlZu\^ = {Z^MIZh, 



al 



a 



(8) 



The remaining 7x7 matrix f [Zj^MpZ h) ^a ) (a, /3 = 2, . . . , 8) can be further 

diagonalized, and then gives seven diagonal masses. The charged scalars 
mass matrix M'^j^ also contains the massless unphysical Goldstone bosons 
G^, which can be written as 



G 



'± 



1 



with the unitary matrix Zh and 
{Zj^Ml^Zu),^ 



:{vdH^ - VuH^ + v^^e%^ 



(9) 



zJ,mi,Zh)^^ 



0, 



a 



(10) 



In the physical (unitary) gauge, the Goldstone bosons G^ and G^ are eaten 
by Z-boson and I^-boson, respectively, and disappear from the Lagrangian. 
Then the mass squared of charged and neutral gauge boson are 
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and 



tan/3 



Vu 



\lv\Vv~v^^ 



Here e is the electromagnetic coupling constant, s^ 
with Q^ is the Weinberg angle, respectively. 



sin Q^ and c^ 



(12) 
cos ^„, 



3. Lepton-flavor violation in the /xi/SSM 

In this section, we present the analysis on the decay width of the rare 
LFV processes /J — )■ l^^ and /J — ;■ l^l^l^ in the /iz/SSM. For this study we 
will use the indices /3, C = 1, • • • , 5, a, p = 1, . . . , 8, and r], a = 1, ... ,10. And 
the summation convention is implied on the repeated indices. 

3.1. Rare decay I J — )■ l^'j 



l-iip) A h(p + q) hip) S„,P„ h{p + q) 
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Figure 1: Feynman diagrams for the LFV process L — ^ Zj 7. (a) represents the contri- 
butions froni neutral fermions x'^ a-nd charged scalars S^ „ loops, while (b) represents the 
contributions from charged fermions xpx ^^'^ neutral scalars Na {N = 5*, P) loops. 

The amplitude for Ij — > l^'j (including /i — >■ 67 and r — )■ ^27 ) is generally 
written as 18 



T = ee'^u^ip + q) [q\iA^PL + A^Pn) 
+ miia^.q^iA^PL + A^Pr)] u,{p) 



(13) 



where q is the injecting photon momentum, p is the injecting lepton momen- 
tum, and mi is the mass of the j-th generation charged lepton, respectively. 
Furthermore, e is the photon polarization vector, Ui{p) {vi{p) in the expres- 
sions below) is the wave function for lepton (antilepton), and P^ = |(1 — 75), 



Pr = |(1 +75)- Here, the Feynman diagrams contributing to the above am- 
phtude are shown in FigUl And the coefficients can be written by 

j^L,R ^ ^{n)L,R ^ j^{c)L,R ^^ = 1,2), (14) 

where Aa ' denote the contributions from the virtual neutral fermion loops, 
and Aa ' stand for the contributions from the virtual charged fermion 
loops, respectively. After integrating the heavy freedoms out, we formulate 
those coefficients as follows 
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(15) 

where the concrete expressions for form factors Ik {k = 1,...,4) can be 
found in Appendix E Additionally, x = m'^/m^^, m is the mass for the 



corresponding particle and mw is the mass for the PF-boson, respectively. 
In a similar way, the corrections from the Feynman diagrams with virtual 
charged fermions are 

1 
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(16) 



Using the amplitude presented in Eg. (1131) . we then obtain the decay width 



for / . — )■ Z^ 7 as |18| 
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And the branching ratio of / • — ?■ /j 7 is 



Br(/- ^ l;i) 



m -> ^r7) 



(17) 



(18) 



where F^- denotes the total decay rate of the lepton Z • . In the numerical 
calculation, T^ ^ 2.996 x lO-^^^GeV for the muon and T^ ^ 2.265 x lO^^^Qgy 
for the tauon. 

3.2. Rare decay Ij — )■ l^l'^lf 
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Figure 2: Penguin-type diagrams for the LFV process l^ — > l^l^li in which a photon 7 
and Z-boson are exchanged. The blob indicates an l^ ^ K ~ ^ vertex such as FigH] or 
l~ ~lj ~ Z vertex where the Z-boson is externaL 



For the rare LFV processes /J -^ KK^t (including \i — )■ 3e), the cor- 
responding effective Hamilton originates from penguin-type diagrams and 
from box-type diagrams. The 7-penguin contribution can be computed us- 
ing Eq. flT3|) . with the result 



T^-^ = u,{vx) q'lMiPL + A'^PR) + mi^ia^,q''{A^PL + A^PR) Uj{p) 



X -^Ui{p2)j^Vi{p3) - (pi ^ P2) 



(19) 



Similarly, the contribution from Z-penguin diagrams which are depicted by 
FigElis 
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where mz is the mass for the Z-boson and 
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The contributions to the effective couphngs -F^^ and Fj^j^ are 
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Here, the concrete expressions for Gk are given in Appendix E 
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Figure 3: Box-type diagrams for the LFV process h — > l.^ h ^t ■ (^) represents the contri- 
butions from neutral fermions Xn o- ^^id charged scalars S~ „ foops, and (b) represents the 
contributions from charged fermions xp,c ^-nd neutral scalars -/Va,p {N = S, P) loops. 



Furthermore, the effective Hamilton from the box-type diagrams which 
are drawn in Figl3]can be written as 
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with 
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(24) 



The effective couphngs Ba ' originate from those box diagrams with virtual 
neutral fermion contributions: 
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Correspondingly, the effective couplings from the box diagrams with virtual 
charged fermion contributions Ba ' are 
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Using the expression for the above amphtude, we can calculate the decay 



width for /" -> /-/r/+ 18 
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And the branching ratio of l^ — )■ l^ l^ if is 
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4. (g - 2)^ in the fiuSSM 

The anomalous magnetic dipole moment (MDM) of the muon can be 
actually be written as the operator 



^MDM — T O-filfiC" l^iFap , 



(30) 



where a"^^ = |[7",7'^], i^a^ is the electromagnetic field strength, l^ denotes 
the muon which is on-shell, m^ is the muon mass and a^ = ^{g — 2)^. 
Adopting the effective Lagrangian approach, we can get [19 
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where Qf = —1, 5R(- ■ ■) represents the operation to take the real part of a 

r D 

complex number and Cg g denote the Wilson coefficients of the corresponding 
operators O2 g 
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In the /iz/SSM, the SUSY corrections can be written as 
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The effective couplings Cg g represent the contributions from the triangle 
diagrams with virtual neutralinos 
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Similarly, the contributions C2 q originating from triangle diagrams with 
virtual charginos are 
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5. The numerical results 

5.1. The parameter space 

It is well known that there are many free parameters in various SUSY 
extensions of the SM. In order to obtain a more transparent numerical re- 
sults, we take some assumptions on parameter space of the /xz/SSM before 
we perform the numerical analysis. 

In lepton sector, we adopt the minimal flavor violation (MFV) assump- 
tions 

i^ijk = 1^ and {A^K)ijk = A^k, ii i = j = k, and zero otherwise. 



m, 



2r 2 2r 2 2r 

m^dij, m-c = m-cdij, m-c = m-cdij, 



Yu^j — Yi/iOij, Ye^j — Ye^Oij, Aj — A, Vuf — Vyc, 

{AMij = A,Y,^5,j, {A,Y,),j = A,YJ,j , and (AaA). = Ax\ (36) 

where i, j, k = 1, 2, 3. 

The 3x3 matrix Yy determines the Dirac masses for the neutrinos Y^jVu ~ 
rriDi and the tiny neutrino masses are obtained through TeV scale seesaw 
mechanism itli, ~ mDm~^nijy. This indicates that the nonzero VEVs of 
left-handed sneutrinos satisfy v^. ^ Vu,d, then 



v„ 



tan/3 ~ — . 

Vd 



(37) 
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Assuming that the charged lepton mass matrix in the flavor basic is in 
the diagonal form, we get 



Y -^ 

Vd 



(38) 



where m/^ is the charged lepton /j mass, and we parameterize the unitary 
matrix which diagonalizes the effective light neutrino mass matrix rrieff (can 
be found in [Appendix Cp as [20 



f/. 



C12C13 

-S12C23 - Ci2S23Sl3e' 
S12S23 - Ci2C23Si3e*'' 

X diag{l,e^ 



^21 ■■ "31 ■ 
2,6 2 



'5l2Cl3 ^136 

C12C23 - Si2S23-Sl3e*'^ S23C13 
-C12S23 - Si2C23Si3e*'^ C23C13 



(39) 



where Cjj 



cos 6ij, Sij = sinOij, the angles 6'^ = [0,7r/2], 5 = [0,27r] is 
the Dirac CP violation phase and 021, 031 are two Majorana CP violation 
phases. Here, we choose 6 = a^i = asi = 0. Ui, diagonalizes riieff in the 
following way: 



U^mlfjUieffU^ 



diag{ml ,ml , m 



"u-i/ 5 



(40) 



where the neutrino mass niy^ connected with experimental measurements 
through 



ml^ - ml^ 



Am2i, 



ml^ - ml^ 



Aml^. 



(41) 



The combination of Eq. (l39|) . Eq. (l40|) . Eq. (HT]) with neutrino oscillation ex- 
perimental data gives some strong constraints on relevant parameter space 
of the /ii/SSM. 

At the EW scale, the soft masses m% , m% , m| and mL are derived 
from the minimization conditions of the tree-level neutral scalar potential. 



which are given in Appendix A Implying the approximate GUT relation 

Ml 



7;:^^ 



(2 ~ O.5M2, the free parameters affect our analysis are 



A, K, tan/3, Ax,K,u,e, rrigc, v^c, M2 . 



(42) 



To obtain the Yukawa couplings Y^. and v^^ from Eq. fl40|) . we assume 
the neutrinos masses satisfying my^<my^<my^, and choose rriy^ = 10~^eV as 
input in our numerical analysis. Then we can get m,^^ 3 from the experimental 
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data on the differences of neutrino mass squared. For Ui,, the values of 6ij 
are obtained from the experimental data in Eq.O- And the effective light 
neutrino mass matrix me// can approximate as [16l | 



"^e//,^ 



2Av^c 

1 — Oibj 

3A ^ 



1 - 3^jj 



CliO/j 7 



(43) 



where 



,2n2 



A = X\vi + <) + AXnvtcVdVu - UX'v^cAB , 

A = Kvlc + XVdVu , 

1 e^ e^ 

+ 



B c^ Ml s2 M2 ' 



W 



W 



at = Yy^Vu , hi = Yy^Vd + ^Xvu^ . 



(44) 



Then, we can numerically derive Y^^ ~ (9(10 '^) and t;,^. ~ 0(10 ^GeV) from 
Eq.dlOD. 



5.£ Branching ratio of LFV processes 



10-" 




Figure 4: Branching ratio for the process /i — !• 67 varies with M2 for tan/3 = 3, 10, 30, 
respectively. 

Considering the research of the /ii/SSM [^, we choose the relevant pa- 
rameters as A = 0.1, K = 0.01, m^c = A^. = Ax = 1 TeV, A^ = A^ = —1 TeV 
and Vyc = 800 GeV in next numerical analysis for convenience. With those 
assumptions on parameter space, we present the branching ratio of /i — )■ 67 
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Figure 5: Branching ratio for the process ^ — >■ 3e varies with M2 for tan/3 — 3, 10, 30, 
respectively. 



versus M2 in FigJH As M2 < 2 TeV, the theoretical evaluations exceed the 
upper experimental bound easily. The fact implies that experimental data 
do not favor small M2. Along with increasing of M2, theoretical evaluation 
on the branching ratio of /i — )■ 67 decreases steeply. As M2 = 3 TeV and 
tan/3 = 10, theoretical evaluation on the branching ratio of /i — )■ 67 is about 
5 X 10^^'^ which can be detected in near future. In the future, the expected 
sensitivity for Br(/i — )• 67) would be of order 10~^^ J2l| . Differing from LFV 
processes which are researched in the BRpV model J22|, the large VEVs of 
right-handed sneutrinos in the /ii^SSM induce new sources for lepton-flavor 
violation. So, here the branching ratio of ;U — )■ 67 can easily reach the upper 
experimental bound 2.4 x 10~^^ |13J |. 

We also investigate the // — )■ 3e processes in detail. And the branching 
ratio of yU — )■ 3e is also decreases with increasing of M2, and raises with 
increasing of tan /3, which is presented in the FigJS] By Introducing the right- 
handed sneutrinos which the VEVs are nonzero to the yUz/SSM, the branching 
ratio of /i — 7- 3e can also easily reach the upper experimental bound 10~^^ 



13|. We can see that the experimental bounds of the branching ratio of 
/i — )■ 3e and /x — )■ 67 give very strong constraints on the /xi^SSM. 

In Figjni we show the branching ratio for r — )■ /i7 versus M2 as tan/3 = 
3, 10, 30. Similar to the case of /i — )■ 67, the evaluation on the branching 
ratio for r — )■ /i7 decreases with increasing of M2, and is enhanced by large 
tan/3. As M2 = 3 TeV and tan/3 = JXI, Br(r — )■ 1^7) ^ 10""^^ is four orders 



below the expected sensitivity 10 ^ [23 
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Figure 6: Branching ratio for the process t -^ fi^ varies with M2 for tan/3 = 3, 10, 30, 
respectively. 



5.3. Muon anomalous magnetic dipole moment 



l.xlO"* 



8. X 10^' - 



6.x 10"' - 




Figure 7: The SUSY contribution to the anomalous magnetic dipole moment of the muon 
varies with M2 for tan/? = 3, 10, 30, respectively. The gray area denotes the Aa^ at 1.8 
standard deviation. 

Finally, we analyze the anomalous magnetic dipole moment of the muon 

in the /xz/SSM. Rescaled the final result of the E821 Collaboration at BNL 

24| using n/p magnetic moment ratio of 3.183345137(85) from ref. [25|], the 

PDG Collaboration [13.] gives the world average of muon anomalous magnetic 

dipole moment 



cxp 



^{g. 



2) = 11659208.9(5.4)(3.3) x 10"^° 
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(45) 



where the statistical and systematic uncertainties are given, respectively. 
And the Standard Model (SM) prediction [l3| is 



af^ = 11659184.1(4.8) x 10"^° . (46) 

So, the difference between experiment and the SM prediction 

Aa^ = a;fP - af^ = 24.8(8.7) (4.8) x 10"^° , (47) 

represents an interesting but not yet conclusive discrepancy of 1.8 standard 
deviation. An alternate interpretation is that Aa^ may be a new physics sig- 
nal with supersymmetric particle loops as the leading candidate explanation. 
If treated the supersymmetry as the leading explanation, parameter space of 
the /iz/SSM should be constrained by the experimental data on Aa^. 

The SUSY contribution to the muon anomalous magnetic dipole moment 
in the /iz/SSM is shown in FigJTl The result shows that when tan/3 = 3, Aa^ 
constrains M2 < 1 TeV, which is opposite to what the upper experimental 
bound of Br(yU — )■ e7) constrains. The fact implies that experimental data do 
not favor small tan/3 in the ;Uz/SSM with the MFV assumptions ( 136|) . When 
tan/3 = 30, Aa^ constrains 2 TeV < M2 < 7 TeV, compared with that the 
upper experimental bound of Br(/i — )■ e7) constrains M2 > 3.5 TeV, the M2 
has more consistent interval. So, under the MFV assumptions, the /iz^SSM 
favors large tan /3 and M2 for consistent with experimental data. 

6. Conclusions 

Besides the superfields of the MSSM, the /xz/SSM introduces three exotic 
right-handed sneutrinos i)'^ to solve the fi problem of the MSSM. And exotic 
right-handed sneutrinos which the vacuum expectation values are nonzero 
induce new sources for lepton-flavor violation. In addition, from the scalars 
for the /iz/SSM we strictly separate the Goldstone bosons, which disappear 
in the physical gauge. 

Considering the updated experimental data on neutrino oscillations, we 
analyze various LFV processes and {g — 2)^ in the /xi^SSM. Numerical results 
indicate that the new physics corrections dominate the evaluation on the 
branching ratios of LFV processes in some parameter space of the /iz/SSM. 
And the theoretical predictions on the branching ratios of LFV processes 
/i — )■ 67 and /i — )■ 3e for large tan /3 can easily reach the present experimental 
upper bounds and be detected in near future. Additionally, the present 



experimental observations on {g — 2)^ also give very strong constraint on the 
model. Under the MFV assumptions (136 p . the /zi/SSM favors large tan/3 and 
M2 for consistent with experimental data. Certainly, a neutral Higgs with 
mass nih,, ~ 124 - 126 GeV reported by ATLAS H and CMS ^] also 
contributes a strict constraint on relevant parameter space, we will discuss 
this problem elsewhere. 
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Appendix A. Minimization of the potential 

First, the eight minimization conditions of the tree- level neutral scalar 
potential are given below: 



m 



/^2 
2 „, I ^ {„,2 „,2 



+ {XiXjVi^cv^c + XiXivl)vd - Y^^^VuXXkVu^v^c + A^i;^) = , (A.l) 



q2 

+ {XiXjV^^v^^^ + XiXiV^^Vu + Yy^^VyXi^ijkV^^v^^^ - 2XjVdVu) 

- XjKijkVdVu-Vu- + {Yu^Yu^.v^fV^c + Y^^^Y^^^Vu,v^^)vu = , (A.2) 

ry2 

^%^^3 + ^(^rf -^1+ Vu,V^,)Vu, + {A^Y^)..VuV^c + Y^^^KijkVuVu-V^c 

- Y^^^XkVycVycVd - Yy^^Xjvlvd + Y^^Y^^^Vy^V^cVul 

+ Y.^J.^.vlv,^ = , (A.3) 

mlcv^c + {A^Yy) ..Vy^Vu - {AxX)iVdVu + {Ani^)i^,,Vy<^Vyi - 2XjKijkVdVuVyi 

+ XiXjVycivl + vl) + 2KiimKijkVyc^VycVyc + 2Yy^^^KiklVuVy^Vyc 

- Yy^^XkVy^VycVd - Yy^^XiVy^VycVd + Yy^Yy^^Vy^Vy^Vyc 

+ Y^.yu,,vlvy. = , (A.4) 

where G'^ = gj + gl and g^c^ = g2S^ = e. 

19 



Appendix B. Mass Matrices 

In this appendix, we give tlie mass matrices in tlie /xz/SSM. 

Appendix B.l. Scalar mass matrices 

For tliis subsection, we use the indices i,j,k,l,m = 1,2,3 and a = 
1,...,8. 

Appendix B.1.1. CP-even neutral scalars 

In the unrotated basis S"^ = {hd,hu,{i'i)^,{i'i)^), one can obtain the 
quadratic potential 

y quadratic ~ 7^^ MgS . (B.l 

And the expression for the independent coefficients of Mg are given in detail 
below: 

q2 

Ml^h^ = rn]j^ + -^{"ivl -vl + v^^v^J + KXjV^^Vu- + KKvl , (B.2 

fi2 

- 2Y^^^\jVdV^^ + Y^^Y^^^v^cv^c + Y^^^Y^^^v^^v^^ , (B.3 

ry2 

^h^hu = -{Ax\)iVuf - —VdVu + 2\i\iVdVu - \K.ijkVy-^Vyc 



■■-I] 
y2 



2Y,^^XjVuV,^ , (B.4 

a 

2 

2 



^Li>>^)^ = -l^^dVu, - Y^^^{Xjvl + XkV^cV^c) , (B.5 

q2 

^L(u.)^ = -^^uVu, + {A^Y^)^.v^. - 2Y^^\pdVu + Y^^^Kij^v^^v^c 

^lY^^Y^^^v^v,^, (B.6 

^hi{yl)R = -{AxX)iVu + 2XiXjVdV^c - 2XkKijkVuVyc 

-(K^,Afc + y;^,,A,Ki;,e, (b.t; 

^L(i>f)« = -{AxX)iVd + {AyYy).^Vy^ + 2XiXjVuVi,c - 2XkKijkVdVu<r 

+ 2Yy^^KiikVy^Vyc + 2Y^^^Y^^^VuVu<^^ , (B.8 

fi2 fi2 
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M? 



+ (n.^.y;,, + y;,n,jt;.,i;.e , (b.io) 

- (>^.,A, + y;,,A,)t;dt;., + n«(>^,,t^^ + n,,t^.,t^.j . (B.ii) 

We can use an 8 x 8 unitary matrix Rs to diagonalize the mass matrix 

s 

RlMlRs = {Mf^'f . (B.12) 
By unitary matrix Rs, S'^ can be rotated to the mass eigenvectors Sa- 

hd = Rs^Sa, hu = R^Sa, {Vi) = Rs ^a, {^i) = Rs Sa ■ (B.13) 

Appendix B.1.2. CP-odd neutral scalars 

In the unrotated basis P'^ = {Pd,Pu, {hY , {^iY)) oi^^ can also give the 
quadratic potential 

%w„,„ = ip-A4P'. (B.14) 
and the concrete expression for the independent coefficients of Mp 

^PaPd = "^L + ^(^d - ^n + ■^u,V•y^) + XiXjVufV^- + XiXivl , (B.15) 

ry2 
^PuPu = "^k - ^(^d - ^n + Vv^^vr) + \\V^-V^c + Wvl 

- 2Y^,,\jVdVu, + Y^^Y^^^v^cv^c + Y^^^Y^^^v^^v^^ , (B.16) 

^PdPu = iAxX)iV^c + XkKijkV^cv^c , (B.17) 

^hiu.y = -Y..Y>^ivl + XkV^^v,.) , (B.18) 

^PuihY = -i'^Mij'^i^'', - yu^ki^ijkVufV^c , (B.19) 

Mp^(-c)7 = {Ax\)iVu - 2\kKijkVuVu^ - (Yi,^,h - Yi,^i^\i)v^^v^c , (B.20) 
^Pu(i>fy = iAx\)iVd - {AyYy).^Vy^ - 2{\kKukVd - Yy^^KiikVy^)vy<^ , (B.21) 
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+ yu^J.,^Vuiv,^ , (B.22) 

- {Y^^Y,^^ - y;,n,Jt;,,i;,c , (B.23) 

- {Yu^.K + Y,^^\,)vdv,, + Y,^XYu,,vl + nyt;.,i;.J . (B.24) 
Using an 8 X 8 unitary matrix Rp to diagonalize the mass matrix M"^ 

RlMJ,Rp = {Mf'^'f , (B.25) 
we can obtain the mass eigenvectors Pq,: 

Pd = RpPai Pu = RpPa, {i^i) = Rp Pa, {^^i) = Rp Pa ■ (B.26) 

Appendix B.1.3. Charged scalars 
The quadratic potential includes 

Vquadratic = S Mg±S , (B.27) 

where S'^'^ = {H^, H^,ef_,ep_) is in the unrotated basis, e^, = Cj and 
e^ = ef. The concrete expression for the independent coefficients of M|± 
are given below: 



„2 (^2 

2 1^2/2 \ , ^ /- 2 2 



d d ^ 4J: 

+ Fe,,i;^,i^.,i^.,, (B.28) 

„2 ^2 

+ n^,y.,^i;.ct;,., (B.29) 

+ y;,,Ajt;„i;,, , (B.30) 

't% 2 

i 

2 



^J±e± = y^dt^i^. - Yu,^\kV^<i^V^<. - Ye^Y^^^VdVu^ , (B.31) 

^2 
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^Uik^Ukj'^u'^Uj ) 



Ml- 






Mi 



Mi 



1 



^i^ 



{AeYe)ijVd - Ye^^XkVuVyC , 
1 

2 



^k e± = ^%. - ^aliv^ -vl + Vy^VyJSij + Y^^Ye^y^ 



'^e.u^ekj'^i^k'^i^l ■ 



(B.32) 
(B.33) 

(B.34) 



(B.35) 
(B.36) 



(B.37) 



M2 



Through an 8 x 8 unitary matrix Rs± to diagonahze the mass matrix 



S'^ can be rotated to the mass eigenvectors S*^: 



(B.38) 



m 



■)2a c± ~± 



(2+J)a r,± ~± 



,(5+J)a Q± 



JDi-a Q± TT± _ T32a q± 2± — 0^^+*^" C± p;± — Dl.o+*i" C± (Tl. 'i.Q\ 



Appendix B.2. Neutral fermion mass matrix 

Neutrinos mix with the neutrahnos and therefore in the unrotated basis 



loT 



X 



terms in the Lagrangian: 



where 



B°,W°, HdyHu,^^^.,^^.], one can have the neutral fermion mass 



/oT 



M„ 



x'°^M„x'° + H.c. 



M m^ 

m 03x3 



(B.40) 



(B.41) 



with 



Sl 



92 



m 



91 



92 



I 



.£Ln, 92_ Y V c Y V Y V Y V 



■ l^li 

/ 

/ 



^'^'^■i V^^^a U -'i^Si^^'f iv'j.l'^U i^U-j,2''^U J'^l^ss'^M 



■ 1^21' 

■ 1^31 "* 



1^23 '-'« 



(B.42) 
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and 



M 



( Ml 





M2 





91 „, 
-92 „, 





-Ait;„ 





-A2t^« 





-Ast^n 


91 ,, 


-92 „, 


— Ajfjyc 





1/1 


Z/2 


1/3 








-Aii^„ 


yi 


2KujV,yC 


2Ki2jVi,c 


S/tlSj"^;/ 








-A2t^M 


y2 


2K21jV^c 


2K22jVyC 


2k,23jVu^ 


I 





-hvu 


1/3 


2 K,31 jVu^ 


2k,32jVu^ 


2k,33jVu 



(B.43) 



where |/j = —XiVa + Fjy^ti,y.. Here, the submatrix m is neutrahno-neutrino 
mixing, and the submatrix M is neutrahno mass matrix. This 10 x 10 sym- 
metric matrix M„ can be diagonahzed by a 10 x 10 unitary matrix Z„: 



ZlMr,Z„ 



M, 



nd 



(B.44) 



where Mnd is the diagonal neutral fermion mass matrix. Then, we have the 
neutral fermion mass eigenstates: 



with 



/(-a 1 ,^o 



a 



1,...,10 



(B.45) 







Hd 









7{7+«)Q!^o 



(B.46) 



Appendix B.3. Charged fermion mass matrix 

Charged leptons mix with the charginos and therefore in the unrotated 
basis where \E'~^ = ( — iA^, /J^, e^. j and ^+"^ = ( — iA+, iJ^, Cp, I, one can 
obtain the charged fermion mass terms in the Lagrangian: 



^"'Mc^+ + H.c. 



(B.47) 



where 



Mr 



M± h 

C ITLi 



(B.48) 
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Here, the submatrix M± is chargino mass matrix 

M^=( ^' f ^« "l . (B.49) 

And the submatrices h and c give rise to chargino-charged lepton mixing. 
They are defined as 





'i-ei^Vvi ~^ei2'^Ui '~'^eiz'^Pi 




(B.50) 



(B.51) 



And the submatrix irti is the charged lepton mass matrix 

Ye^^Vd Ye^^Vd Ye^sVd • (B.52) 

^e-siVd Ye,^Vd Ye,,Vd J 

This 5x5 mass matrix Mc can be diagonahzed by the 5x5 unitary matrices 
Z_ and Z+: 

Z^M,Z+ = M,d , (B.53) 

where Med is the diagonal charged fermion mass matrix. Then, one can 
obtain the charged fermion mass eigenstates: 

Xa=(^), a = l,...,5 (B.54) 



with 






(B.55) 
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Appendix C. Approximate diagonalization of mass matrices 

Appendix C.l. Neutral fermion mass matrix 



If the R-parity breaking parameters are small in the sense that for 16|, |28 



i = m.M-^, (C.l) 

all ^ij <^ 1, one can find an approximate diagonalization of neutral fermion 
mass matrix. In leading order in ^, the rotation matrix Z„ is given by 



1 - Wi -e \(v 

^ 1 - w )\^u, 



IccT 1 ( n u ) • (C.2) 



The first matrix in (lC.2p above approximately block-diagonalizes the matrix 
Mn to the form diag (M, nieff), where 



m, 



eff 



-m.M-\m^ . (C.3) 



The submatrices V and Ui, respectively diagonalize M and me// in the fol- 
lowing way: 



V^MV = Md , 
UjnieffUu = niyd 



(C.4) 



where M^ and m^d are respectively diagonal neutralino and neutrino mass 
matrix. 

Appendix C.2. Charged fermion mass matrix 

Similarly to the approximate diagonalization of the neutral fermion mass 
matrix discussed above, it's also possible to find an approximate diagonal- 
ization procedure of the charged fermion mass matrix for the small R-parity 
breaking parameters 28|. Then, we can define 






± ' (C.5) 



All $,L,j ^ 1 and ^r.^ <^ 1, so in leading order in ^l and ^r, the rotation 
matrices Z_ and Z+ are respectively given by 

-(-l5S--&)(^-;j. <«) 
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IfT 






2^R^R 
^R 1 ^ 2^R^R 











(C.7) 



Then the matrix M^. can approximately be block-diagonahzed to the form 
diag {M±,mi). And the submatrices U^,U+ and V^-,V"+ respectively diago- 
nalize M± and mi in the following way: 



V^miV+ 



M. 



±d 



mid, 



(C.8) 



where M±d and m;^ are respectively diagonal chargino and charged lepton 
mass matrix. 



Appendix D. Interaction Lagrangian 

In this part, we give the interaction Lagrangian of the relative vertices 
for the LFV processes in the /iz/SSM. And we use the indices i, j = 1, . . . , 3, 
/3, C = 1, • • • , 5, a, p = 1, . . . , 8 and r^ = 1, . . . , 10. 

Appendix D.l. Charged fermion-neutral fermion-gauge boson 

We now give the interaction Lagrangian of charged fermion, neutral fermion 
and gauge boson, 



^X(Xp^ 



yZx(Xp^ 






+ w-ucY'^'^'^rPL 



c 



WX%Xl3 n 



R 



i'Pr)x, + ---, (D-1) 



where the coefficients are 



C 



ZxcXp 



C 



C 



ZxcXp 
R 

WxpX°r, 



'^^W^W 



^^w^w 



;i _ 2s2 )<5«/^ + Z^h Z 



K* yW 



2C* ^2/3 



2zTzl^ + zTzf - 2sU^^ 



+ 



^w- 



c 



c 



Wxpx°r, 

R 

Wx°r,Xp 



V2s, 
e 

711 



\/2Z'Jzr + Z'JZI''* + zi2+^)/^47+.).* 



V2Z'fzl^ - Z'fZ^^ 



c 



WxpX^n' 



c 



Wx°r,Xp 

R 



c 



Wx(iX°n 
R 



(D.2) 
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Appendix D.2. Charged scalars-gauge boson 

The interaction Lagrangian of charged scalars and gauge boson is written 
as 



A„t = teF,S-*d^S- + ieC^''-^^*Z,S:*d^S- + 
The coefficient is 



/^ZSa Sp 



■^^w^w 



(5+i)"* d(5+«)p 



'l-2s'J5''P-R'^r' R's± 



(D.3) 



(D.4) 



Appendix D.3. Charged fermion-neutral fermion-scalars 

The interaction Lagrangian of charged fermion, neutral fermion and scalars 
is similarly written by 



^SaXpXC 



— * — 0/'/^ "^ X/3X77 7") I y^-^a \pXrj T-) \ I 



And the coefficients are 



(D.5) 



C 



SaXpXi 



— e 



V2s 



Rl'^Z'Jzl^ + R]?Zl^Zi^ + i?(f+^)"zi'+^)^z5' 



p(5+J)a^l/3^(2+j)C _ pla^{2+J)/3^(2+i)C 



1 V p(2+j)a 7(2+J)/3 72c ^ x p{2+J)a ^2,3 72C 



c 



-PaX/JXC 



«e 



V2s, 



_^2q^1/3^2C _^ pla^2/3^1C _^ p(5+J)a^(2+i)/3^1C 



+ V^^"^' 



p(5+J)a 7I/3 y{2+3)C „la 7{2+J),3 7(2+i)C 

-ftp ■^_ -^-f — JTp Z;__ ■^ + 



,(2+j)a 7(2+J)/3 ^2C ^ ^ p(2+J)a ^2/3 72C 



- ^y..A 



Zj _ ^ -\- — — "^^i^p Zj _ Zj \ , 

v2 



C 



Sa X°Xl3 



y2swCw 



rj2a* /72/3 



ch/Z^'' + swZ, 



Irj 



Sw 



r>2a* y^P y^V 



L^i 



Sw 



p(5+.)a*^(2+.)/3^1, ^ y^^^J?(2+^)"^2/3^(4+,), 



28 



(2+i)/3 



C 



Sa*Xl3X?i 






y/2swCw 



■^^a ^(7+.), _ i?(2+0"^3J _ X,RllZfzi^+'^\ 



Rlrzl^ + R 



{2+i)a* ry{2+i)l3' 



Cw 



Zl' + swZl^ 






R'sTzl"^ + i?^'±+')"*z(^+^)''j + y;^^j?|"±zi'+')''z(^+^> 



+ r,^.<±^^'^" 



c 



c 



SaXpXC 

R 

SaX°r,Xli 

R 



C 



SaXCX/S 



^2/3^(7+.), _ ^{2+.)/3^3J _ A,i?|"±Z!^Z(4+*)^ 

^PaXCX~^* 



^PaXfiXC 



'R 



c 



•Sa'X/sxSJ 



C 



Sa*XfiX°r, 



R 



C 



Sa X°r,Xfi 



Appendix E. Loop-momentum integral 

2 

Defining Xi = ^-, we can find tlie loop- momentum integral for /~ 



Il(xi,X2) 
l2{Xi,X2) 
l3{xi,X2) 



l + lna;2 Xi In aji — X2 In 0:2 



IGtt^ \-{x2 — Xi) 



{X2 -XiY 



167r2 
1 



1 + lnxi Xi In Xi — 0:2 In X2 



{X2 - Xi) 



{X2 -XiY 



(D.6) 

(E.l) 
(E.2) 



3 + 21na;2 2x2 + 4:X2\nx2 2x^lnxi 



327r^ L (x2 — Xi) 
+ 



{X2 -XiY 



{X2 - XiY 



2x1 1'^ ^2 



(X2 -XiY 



(E.3) 



/4(a;i,a;2) 



1 



967r2 



ll + 61na;2 15x2 + 18x2 In 0:2 6x2 + 18x2 In X2 



(X2 -Xi) 



(X2 - Xi' 



+ 



(X2 - Xi)^ 



+ 



6xf Inxi — 6x2 lnx2 



(X2 - Xi) 

And we also can find the loop- momentum integral for /~ — )■ l^l^l^: 



(E.4) 



Gi(xi,X2,X3) 



1 



Xi Inxi 



+ 



X2 In X2 



IGtt^ L(xi - X2)(xi - X3) (x2 - a;i)(x2 - X3) 
X3 In X3 

(X3 -Xi)(x3 -X2)J' 



+ 



(E.5) 



G'2(Xi,X2,X3) 



167r2 



-(A + l + lnx^) 
29 



xi Inxi 



(Xi -X2)(Xi -X3) 



+ 



Xo In X2 



x\ In X3 



(X2 -Xi)(x2 -X3) (Xs -Xi)(x3 -X2) 



(E.6) 



2 

Here, x„ = -^. G2ixi^X2-,x-i) is divergence, so here we use dimensional 

regularization to cancel the divergent part (A + 1 + Inx^). In the numerical 
calculation, we will keep the remaining convergent part. 



Gz\Xi,X2,X^,Xi) 



167r2 



Xi Inxi 



(Xi -X2)(xi -X3)(xi -Xi) 



X2 In X2 



+ 



X'i In X3 



(a;2 - a;i)(x2 - 0:3) (0:2 - X4) (xg - Xi)(x3 - X2)(x3 - X4) 
X4 In X4 



(X4 - Xi)(x4 - X2)(X4 - X3) 



(E.7) 



Gii[xi, X2, X3, X4) 



167r2 



x\ In Xi 



(Xi - X2)(Xi - X3)(Xi - X4) 



Xn In X2 



+ 



Xo In X3 



X2 - Xi)(x2 - X3)(X2 - X4) (X3 - Xi)(x3 - X2)(x3 - X4) 



+ 



X4 In X4 



(X4 - Xi)(x4 - X2)(X4 - X3) 



(E.8) 
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